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Abstract. The thermodynamical concept of entropy was introduced by Clausius in 1865 in order to con- 
struct the exact differential dS = 8Q/T, where SQ is the heat transfer and the absolute temperature T 
its integrating factor. A few years later, in the period 1872-1877, it was shown by Boltzmann that this 
quantity can be expressed in terms of the probabilities associated with the microscopic configurations of 
the system. We refer to this fundamental connection as the Boltzmann-Gibbs (BG) entropy, namely (in 
its discrete form) Sbg = — k ^2 i=1 Pi lnpi, where k is the Boltzmann constant, and {pi} the probabilities 

corresponding to the W microscopic configurations (hence Yli=iPi ~ !)■ This entropic form, further dis- 
cussed by Gibbs, von Neumann and Shannon, and constituting the basis of the celebrated BG statistical 
mechanics, is additive. Indeed, if we consider a system composed by any two probabilistically independent 
subsystems A and B (i.e., pf+ B = pfpf, V(i, j)), we verify that S B a(A + B) = Sbg(A) + Sbg{B). If a 
system is constituted by TV equal elements which are are either independent or quasi-independent (i.e., not 
too strongly correlated, in some specific nonlocal sense), this additivity guarantees Sbg to be extensive 
in the thermodynamical sense, i.e., that Sbg(N) oc TV in the TV >> 1 limit. If, on the contrary, the cor- 
relations between the TV elements are strong enough, then the extensivity of Sbg is lost, being therefore 
incompatible with classical thermodynamics. In such a case, the many and precious relations described in 
textbooks of thermodynamics become invalid. Along a line which will be shown to overcome this difficulty, 
and which consistently enables the generalization of BG statistical mechanics, it was proposed in 1988 the 
entropy S q = k[l — X/i=iP?]/(9 — 1) (l $i = Sbg)- In the context of cybernetics and information 

theory, this and similar forms have in fact been repeatedly introduced before 1988. The entropic form 
S q is, for any 9 / 1, nonadditive. Indeed, for two probabilistically independent subsystems, it satisfies 
S q (A + B)/k = [S q {A)/k] + [S q (B)/k] + (1 - q)[S q (A)/k][S q (B)/k] / S q (A)/k + S q {B)/k. This form will 
turn out to be extensive for an important class of nonlocal correlations, if q is set equal to a special value 
different from unity, noted q en t (where ent stands for entropy). In other words, for such systems, we verify 
that Sq ent (N) oc TV (TV >> 1), thus legitimating the use of the classical thermodynamical relations. Stan- 
dard systems, for which Sbg is extensive, obviously correspond to q m t = 1- Quite complex systems exist in 
the sense that, for them, no value of q exists such that S q is extensive. Such systems are out of the present 
scope: they might need forms of entropy different from S q , or perhaps — more plainly — they are just 
not susceptible at all for some sort of thermostatistical approach. Consistently with the results associated 
with S q , the ^-generalizations of the Central Limit Theorem and of its extended Levy-Gnedenko form 
have been achieved. These recent theorems could of course be the cause of the ubiquity of g-exponentials, 
q-Gaussians and related mathematical forms in natural, artificial and social systems. All of the above, as 
well as presently available experimental, observational and computational confirmations — in high energy 
physics and elsewhere — , are briefly reviewed. Finally, we address a confusion which is quite common in 
the literature, namely referring to distinct physical mechanisms versus distinct regimes of a single physical 
mechanism. 

PACS. 05.20.-y Classical statistical mechanics; 02.50.Cw Probability theory; 05.90.+m Other topics in 
statistical physics, thermodynamics, and nonlinear dynamical systems; 05. 70. -a Thermodynamics 

1 Introduction 

* To appear in Statistical Power-Law Tails in High Energy 
Phenomena, ed. T.S. Biro, Eur. Phys. J. A (2009). The concept of entropy S, as well as its name, were intro- 

a tsallis@cbpf.br duced in thermodynamics by Clausius in 1865 pQ. It was 
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done on a purely macroscopic basis (the possible existence 
of a microscopic world was merely speculative at the time, 
with just a few incipient scientific indications), with the 
purpose to associate with the heat transfer SQ an exact 
differential. This was indeed achieved through the cele- 
brated relation dS = 5Q/T, where dS is the differential 
entropy, and the absolute temperature T the integrating 
factor. 

A decade later, in the period 1872-1877, it was shown 
by Boltzmann [2j[3] that this quantity can be expressed 
in terms of the probabilities associated with the micro- 
scopic configurations of the system. We refer to this con- 
nection — one of the deepest ever done in physics — as the 
Boltzmann-Gibbs (BG) entropy. In its present (discrete) 
form, it is written as follows 



Srg = ~ky i Pi In pi 



(1) 



where k is the Boltzmann constant (or some other conve- 
nient value, e.g. k = 1, in areas outside physics, such as 
information theory, cybernetics and others), and the prob- 
abilities {pi} corresponding to the W microscopic config- 
urations satisfy 

w 

X> = 1 - (2) 

i=l 

This entropic form, further discussed by Gibbs [1], von 
Neumann [5] and Shannon [BJ, constitutes the basis of 
the BG statistical mechanics, one of the monuments of 
contemporary physics. Eq. (JTJ satisfies a variety of conve- 
nient mathematical properties (non-negativity, concavity, 
expansibility, Lesche-stability, composability, Topsoe fac- 
torizability, finite entropy production per unit time satis- 
fying the Pesin identity). For equal probabilities, i.e., for 
Pi = l/W(Vz), it attains its maximal value, namely 



Sbg = k\nW . 



(3) 



In the present context, let us focus on its additivity prop- 
erty. An entropy S({pi}) is said additive [7] if, for a system 
composed by any two probabilistically independent subsys- 
tems A and B (i.e., satisfying p. 



A+B 



pfpf, V(i,i); i 



1,2,..., W A ; j = 1,2,..., W B ), we verify 



S(A + B) = S{A) + S(B), 



(4) 



where S(A + B) = S{{pf+ B }), S{A) = S({pf}) and 
S(B) = S({p?}). It is straightforward to verify that 
Sbg({k}) given by Eq. (fTJ) is additive. Due to this prop- 
erty, the BG entropy of any system made of N equal and 
independent elements satisfies 



Sbg(N) = NS BG (l) 



(5) 



This fact obviously complies with the classical thermody- 
namical requirement for the entropy S to be extensive, 
i.e., such that 



,. S(N) 

hm — — - < oo . 



Indeed, in such a case, 

lim = s BG (l) < kln[W(l)\ , (7) 

AT— >oo TV 

where W(l) is the number, assumed finite, of possible con- 
figurations of one element. 

If the system is constituted by N equal elements which 
are not strictly independent, but quasi-independent in- 
stead (i.e., not too strongly correlated, in some nonlo- 
cal sense to be further clarified later on; typically for a 
Hamiltonian many-body system whose elements interact 
through short-range interactions, or which are weakly quan- 
tum -entangled) , the additivity of Sbg guarantees its ex- 
tensivity in the thermodynamical sense, i.e., that Eq. ([6j 
is satisfied. 

If, on the contrary, the correlations between the N ele- 
ments are strong enough (a feature which might typically 
occur for nonergodic states, e.g., in Hamiltonian many- 
body systems with long-range interactions, or which are 
strongly quantum-entangled) , then the extensivity of Sbg 
might be lost (at least at the level of a large subsystem 
of a much larger system), being therefore incompatible 
with classical thermodynamics. In such a case, many of 
the useful relations described in textbooks of thermody- 
namics may become invalid. It is precisely this pathologi- 
cal class of systems the one which is addressed within the 
thermostatistical theory usually referred to as nonexten- 
sive statistical mechanics [5 |ll8pi9| . described in the next 
Section. 



2 Nonadditive entropy and nonextensive 
statistical mechanics 

2.1 Nonadditive entropy S q 

As an attempt to generalize BG statistical mechanics, and 
possibly provide a frame for handling some of the above 
mentioned pathological systems, it was postulated in 1988 
[8] the following entropy: 

1-T W v q 

S = k ^ i=lPl (q e R; Si = S BG ), (8) 

q-l 

A simple manner to obtain Si = Sbg 1S through the use of 
pi' 1 = e (?-i)inw ~ l + {q-l)\n Pl . \{q< 0, the sum must 
be done only over configurations which have nonzero prob- 
ability to occur. Entropy (jSJ) can be conveniently rewritten 
in the following alternative forms: 



w 



(6) 



S Q = k ^2Pi ln q( l /Pi) 

i=l 
W 

= -k^2pfln g pi 

i=l 
W 

= -k^pi ln 2 - q Pi ■ 



(9) 
(10) 

(11) 
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where the q-logarithmic function is defined as follows: 

z l ~ q - 1 

ln„ z = (q E M; z > 0; ln x z = lnz). (12) 

1-9 

S q attains its extremum (maximum for q > 0, and min- 
imum for q < 0) for equal probabilities, and its value is 
given by 



W(N) ~ Cfi N (with C > and (jt > 1), then q ent = 1, 
i.e., Sbg{N) oc N (N » 1). But if strong correla- 
tions forbid many (typically most) microscopic configu- 
rations to occur, then it might happen that the num- 
ber W e ff(N) of effective (or admissible) configurations 
satisfies W e ff(N) « W(N). If we have, in particular, 
W e ff ~ DNP (with D > and p G M), then 



A: In, W. 



(13) 



It can be shown to satisfy, for independent subsystems, 



S q (A + B) 
k 



S q (A) , S q (B) MA^S^)^ 
' '■ k 



k 



Therefore, this entropy is generically nonadditive. It satis- 
fies, nevertheless, most other properties (mentioned previ- 
ously) of the entropy Sbg- In other words, it constitutes a 
sort of minimalistic generalization of Sbg- From Eq. (|14ll 
we obtain 



S q {A + B) = S q {A) + S q (B) 



S q (A)S q (B), (15) 



which exhibits the equivalence between (q — 1) — > and 
k — > oo. Since for stationary states (e.g., thermal equi- 
librium), k appears multiplicatively accompanied by the 
temperature T (i.e., in the form kT), k — > oo turns out to 
be equivalent to T — > oo. We may consider that it is here 
where the fact emerges that, at the T — > oo limit, all the 
microcanonical, canonical, grand-canonical ensembles of 
classical or quantum (Fermi-Dirac, Bose-Einstein, Gentile 
parastatistics) statistics, g-statistics (as will become evi- 
dent later on)), coincide, and coincide with the hypothesis 
of equal probabilities for an isolated system. 

Eq. (fT"4|) can be generalized in the presence of arbitrary 
correlations between two sysbsystems A and B of a given 
system. It becomes [9] 



S q (A + B) _ S q (A) S q {B\A) 



k 



, (1 q) S q (A) S q (B\A) 



S q (A\B) , S q (B) 



, (1 q) S q (A\B) S q (B) 



where S q (A + B) is to be calculated with the joint prob- 
abilities {pfj +B }i Sq(A) with the marginal probabilities 
{P?} = {Y,T= B iP£ +B } (analogously for S q (B)), and 
S q (A\B) with the conditional probabilities {pfj +B /pf} 
(analogously for S q (B\A)): see for full details. Eq. 
straightforwardly recovers Eq. (jT5J) if A and B are inde- 
pendent, hence S q (A\B) = S q {A) and S q {B\A) = S q {B). 

It is precisely this nonadditivity the property which 
enables thermodynamical extensivity. More precisely, if 
both A and B are very large (i.e., Na >> 1 and Nb » 
1), then a value of q, noted q en t, might exist for which 
S Qsnt (A + B) ~ S qmt (A) + S qent (B). In other words, if 
a system has N >> 1 equal elements, it becomes possi- 
ble that a special value of q exists such that generically 
< Iimjy^oo [S qent (N)/N] < oo. 

This interesting feature can be easily illustrated in the 
case of equal probabilities, for which Eq. (fl3|) holds. If 



Qent — 1 • 

P 



(17) 



This type of highly restricted phase space may occur in 
various systems, as has been numerically or analytically 
illustrated in various examples. Let us briefly mention here 
two of them that are analytically tractable, namely an 
abstract probabilistic one and a physical one. 

The probabilistic model consists in N correlated dis- 
tinguishable binary variables [TU] • The probabilities of the 
2 N states vanish excepting for ~ (d+ 1)N of them (which 
can be seen, in the classical Pascal-like triangular repre- 
sentation, as a iV-long "strip" whose width is (d + 1), 
d being non-negative). This model asymptotically satis- 
fies probabilistic scale-invariance (Leibnitz triangle rule) 
in the limit N — > oo. It can be verified that 



Qent 



1 



(18) 



The physical model corresponds to a long ring of N 
1/2 spins with ferromagnetic first-neighbor interactions at 
zero temperature. The interactions are of the anisotropic 
XY ones in the presence of a transverse external magnetic 
field (i.e., along the Z direction) at its critical value. Two 
well known universality classes are contained within such a 
system, namely the Ising universality class (corresponding 
to a central charge c — 1/2), and the isotropic XY uni- 
versality class (corresponding to a central charge c = 1). 
We consider a block of L successive spins among those 
N spins, and address the entropy S q (L) of the N — > oo 
quantum system. More precisely, we are interested in 



(16) 



S q (L) 



l-Tr(p L ) 
Q-l 



(19) 



where p^ = limN^ QO Tr^_Lyp(N) 1 p(N) being the den- 
sity matrix associated with the system of N spins, and 
where we have traced over all but the L successive spins. 
We define, in this case, q en t as the value of q for which the 
block entropy S q (L) is extensive, i.e., such that S q<lnt (L) oc 
L. Such value does exist [11], and it is given by q en t — 
\/37-6 ~ 0.0828 for c = 1/2, and q ent = VlO-3 ~ 0.1623 
for c = 1. By using a recent result within conformal quan- 
tum field theory [T^] , these two values for q ent can be gen- 
eralized for the entire class of (1 + l)-dimensional models 
characterized by a generic central charge c. It is obtained 

El 



Qent 



V9T? - 3 



(20) 



As we see, q en t monotonically increases from zero to unity 
(the BG value!) when c increases from zero to infinity. For 
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Table 1. QSS stands for quasi- stationary state ( [161117] and 
references therein). 

SYSTEM ENTROPY S BG ENTROPY S q (q < 1) 
(additive) (nonadditive) 

Short-range 
interactions, 
weakly 
entangled 

blocks, etc EXTENSIVE NONEXTENSIVE 

Long-range 

interactions 

(QSS), 

strongly 

entangled 

blocks, etc NONEXTENSIVE EXTENSIVE 



c = 4, one obtains q = 1/2, which has already emerged in 
general relativistic problems [13] ; the possible connection, 
if there is one, remains however without explanation at 
the present time. Also, for c = 26, which corresponds to 
string theory [14], we obtain q ent ~ 0.8913. Finally, it 
also remains presently without explanation the reason for 
which c — * oo leads to the BG result, i.e., q ent — 1. 

The previous system is a fermionic d = 1 one, for 
which it is known that Sbg(L) tx InL. Results, though 
only numerical, also exist for a bosonic d = 2 system [11] . 
Once again extensivity only occurs for q ent < 1, whereas 
Sbg(L) oc L. Finally it is known that, for the black hole, 
the BG entropy is proportional to the area L 2 , instead 
of being proportional to the volume L 3 (see, for instance, 
[15]): in fact, more generally, it is known, for (i-dimensional 
bosonic systems, the so-called area law, i.e., the fact that 
Sbg{L) oc L d_1 , which obviously violates classical ther- 
modynamics. All these anomalies are believed to be a con- 
sequence of strongly nonlocal quantum entanglement. In 
fact all of the above results can be unified through the 
following conjectural expression, Vd: 

S BG (L)(xln 2 _ d L= L d 1 _~ 1 ¥=L d cxN. (21) 

In all these cases, as for the d = 1 and the d — 2 above 
described examples, it might well exist a value of q ent < 
1 such that thermodynamic extensivity is ensured, i.e., 
Sq crlt (L) ex L d oc N. This would of course mean that 
Clausius entropy should be, for this class of anomalous 
(sub)systems, identified with S qcnt and not with Sbg- The 
plausible conjectural scenario is summarized in Table 1. 

2.2 Nonextensive statistical mechanics 

Since S q generalizes S bg and maintains most of its math- 
ematically convenient properties (e.g., concavity, Lesche 
stability, among others), it is quite natural to attempt the 
^-generalization of BG statistical mechanics itself. This 
extended theory is usually referred to in the literature as 
nonextensive statistical mechanics. It was first proposed 



in 1988 [8], and later on connected to thermodynamics 
18, 15]. It has since then received a considerable amount 
of applications and verifications in natural, artificial and 
social systems [20 p 21 | l2l ^ l2"3 l 2"4]. Some of its predictions 
have been experimentally and observationally checked in 
systems such as the motion of Hydra viridissima [25] and 
cells [26 ] I27 | I28 | . defect turbulence [29], solar wind [30 1 131 ] 
32 , 33] , cold atoms in optical dissipative lattices pH] , dusty 
plasma [35], silo drainage [36ll37j . high-energy physics (see, 
e.g., [38,39,40,41,42 ). They have also been checked an- 
alytically and computationally in various nonlinear dy- 
namical problems such as the edge of chaos of simple 
unimodal dissipative maps [43] j 44l 45 . 46, 47 ll48H 49| J 50| J 5 H 
52,53 , and long-range-interacting many-body Hamilto- 
nian systems 54,55,16,17 ; also, various applications to 
the so-called scale-free networks are available 56 , 57, 58 , 
59, 60, 61, 62, 63, 64J. In one way or another, most if not all 
of these systems appear to share slow (power-law rather 
than exponential) sensitivity to the initial conditions. In 
other words, for classical systems, at the level of first prin- 
ciples, BG statistical mechanical concepts are legitimate 
and fruitful when the system exhibits a positive maximal 
Lyapunov exponent (corresponding essentially to Boltz- 
mann's molecular chaos hypothesis), whereas vanishing 
maximal Lyapunov exponent appears to be necessary (al- 
though probably not sufficient) for the applicability of the 
nonextensive statistical mechanical concepts. The mecha- 
nisms that typically yield q-statistics involve, at the meso- 
scopic level, non-Markovian processes [65], multiplicative 
noise [66] . nonlinear Fokker-Planck equations 67,68,69, 
[70] . and similar ones. 

Let us now briefly review, within the present theory, 
two important stationary-state distributions, namely the 
g-generalization of the celebrated BG weight (discrete case), 
and the g-generalization of the Gaussian distribution (con- 
tinuous case). 

To generalize the BG factor for the canonical ensemble 
(i.e., a system in a stationary state due to its contact with 
a "thermostat") we follow [19]. We must extremize S q as 
given by Eq. with the constraints 

w 

X> = 1 ' (22) 

i=i 

and 

w 

J2 E i P i q) = u i> (23) 

i=l 

where the escort distribution {P^} is defined through 

P<«> = _EL- p ^ h ( 24) 

and U q is a finite quantity characterizing the width of 
the energy distribution {pi}; {Ei} are the eigenvalues of 
the system Hamiltonian (with the chosen boundary condi- 
tions). Notice, by the way, that constraint (|22]) can equiv- 
alently be written as Y^=i = 1- 



Constantino Tsallis: Nonadditive entropy: the concept and its use 



5 



Through the entropic optimization procedure, we ob- 
tain straighforwardly 



as well as relations such as 



with 



and 



-MEi-U q ) 



Pi — q 



S=lP 3 



ir 



-0q{Ei — U q ) 



(25) 
(26) 

(27) 



with the q-exponential function (inverse of the previously 
denned ^-logarithm) &\ = [1 + (1 - q)z] l/{1 - q) if 1 + 
(1 — q)z > 0, and zero otherwise (ef = e z ), (3 being the 
Lagrange parameter associated with the constraint (|23|) . 
Eq. 1|25[) makes explicit that the probability distribution 
is, for fixed j3 q , invariant with regard to the arbitrary 
choice of the zero of energies. The stationary state (or 
(meta)cquilibrium) distribution (|25[) can be rewritten as 
follows: 

o H 

-1 



Pi 



Z 'a 



with 



and 



Z' 



w 

£< 

i=i 



P'a 



Pa 



1 + (1 - q)0 q U q 



(28) 



(29) 



(30) 



The form (|28[) is particularly convenient for many appli- 
cations where comparison with experimental or computa- 
tional data is involved. Also, it makes clear that pi asymp- 
totically decays like l/IS^ 9 ~ 1 ' for q > 1, and has a cutoff 
for q < 1, instead of the familiar exponential decay with 
Ei for q = 1. 

The connection to thermodynamics is established in 
what follows. It can be proved that 



1 

T 



dSq 
dU a 



(31) 



with T = l/(k(3). Also we can prove, for the free energy, 

1 



F„ 



where 



U q - TS q 



In, Z q = In, Z q 



P 



In, Z q , 



(3U q 



(32) 



(33) 



This relation takes into account the trivial fact that, in 
contrast with what is usually done in BG statistics, the 
energies {Ei} are here referred to U q in Eq. (|25|) . It can 
also be proved 

d 



U a 



d(3 



In, Z q , 



(34) 



r - T dS q 



dU„ 



d 2 F 

dT or 2 



(35) 



In fact the entire Legendre transformation structure of 
thermodynamics is g-invariant, which is both remarkable 
and welcome. 

As a final remark, let us stress an interesting feature 
concerning q > 1 (power-law decay for pi). Let us assume 
that the energy spectrum has a (quasi-continuous) den- 
sity state g(E). The normalization condition ([22]) implies 
that / dE g(E)p(E) is finite. Since, for q > 1, p(E) decays 
as l/E 1 '^- 1 ^, it follows that g{E)/E 1 ^^ must be inte- 
grable at infinity. This determines the maximal value of q 
which is mathematically admissible in the present theory 
(for instance, if g(E) is constant, then q < 2 must be satis- 
fied). Let us now focus on the other constraint, namely Eq. 
(f23|) . We immediately see that Eg(E)[p(E)] q must also be 
integrable at infinity, i.e., Eg(E)/E q /( q ~ 1 '> must be inte- 
grable, which implies the same limit for q as before! So, 
for instance, if g(E) is a constant, both the normalization 
and the energy constraints are finite for q < 2. For q > 2 
the entire theory becomes mathematically inadmissible. 

Let us address now the continuous case which gener- 
alizes the Gaussian distribution. We want to extremize 



S q — k 



1 - / dx [p{x)] q 



with the constraints 



and 



where 



Jd XP ( X = l. 



dxx 2 P {q \x) 



P (q \x) 



Jdx [p(a;)]9 



a 2 being a fixed (positive) quantity. We obtain 



p(x) 



-{3x 2 



J dxe q 



-fix -1 



(36) 



(37) 



(38) 



(39) 



(40) 



where (3 > can be determined by using constraint 
The entire theory is valid for q < 3, above which both the 
normalization and the q-variance (|38p diverge. For q > 1 
the distribution is defined for all values of x; for q < 1 it 
has a finite support. For q < 5/3 the standard variance 
J dxx 2 p(x) is finite; for q > 5/3 it diverges. See [71] for 
a numerical comparison between variance and g-variance, 
which exhibits the considerable convenience of the latter. 



2.3 q-generalized central limit theorems 

We focus here on the q-generalization of the Central Limit 
Theorem (CLT). Let us remind what the standard CLT 
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states essentially. Consider N equal and independent ran- 
dom variables {xj}, j — 1, 2, 3, N, (in fact, quasi- inde- 
pendent random variables, in some specific sense that is 
out of the present brief review, are also admissible) whose 
mean value and variance are finite. For simplicity we as- 
sume that the associated distributions p(x) (one and the 
same, V?) are symmetric (hence the mean value vanishes). 
We are interested in the sum Xn = J2j=i x j (propor- 
tional to the arithmetic average, which is the quantity 
that is typically measured in all generic experiments), and 
in its correponding distribution P(Xjy). The CLT states 
that, after appropriate centering and scaling (with V~N), 
P(Xjy) converges, in the limit N — > oo, onto a Gaussian 
distribution (whose variance coincides in fact with that of 
p(x)). 

Another important theorem, sometimes referred to as 
the Levy-Gnedenko theorem, exists for the case when p(x) 
has a divergent variance (e.g., the Cauchy-Lorentz dis- 
tribution). The random variables are still assumed inde- 
pendent (or quasi-independent, as previously mentioned). 
If the asymptotic behavior of p(x) decays like a/|a;| 1+a 
(a > and < a < 2), then, in the limit N — > oo, 
the attractor is the Levy distribution L a (x) (also called 
a-stable distribution) which, after appropriate centering 
and scaling (with A 1 /"), decays like p{x) itself, i.e., like 
a/\x\ 1+a . 

These two important attractors (in the space of proba- 
bilities) emerge whenever the N variables are independent 
or quasi-independent. What happens if they are not, i.e. 
if they are strongly correlated? The nature of the possible 
attractor depends on the class of correlations. What hap- 
pens then if this correlation is of the q-independent class? 
The precise definition of q- independence (1-independence 
being just plain independence), is given in 72 (see also 
[73]l74] as well as [75]). It is essentially associated with 
the strong correlations which emerge in the presence of 
probabilistic scale-invariance (Leibnitz triangle rule, in 
the case of N binary random variables; see [76j and refer- 
ences therein) Q. 

If the N variables are q-independent then the attrac- 
tors are q-Gaussians with q > 1 if an appropriately q- 
generalized variance of p{x) is finite [72] (see [77] for the 
q < 1 case), and the so-called (q, a)-stable distributions 
if this variance diverges (and p(x) asymptotically decays 
as a power- law) |75j . For q-independent random variables, 
we may summarize the situation as follows: if the appro- 
priately q-generalized variance is finite, the attractor is 
a (/-Gaussian (Gaussian if the variables are independent, 
and q-Gaussians with q ^ 1 otherwise); if it diverges, 
the attractor is a (q, a)-distribution (Levy distributions 
if the variables are independent, and (q, a)-distribution 
with q ^ 1 otherwise) 0. 

1 The specific q-independent correlation is based on the q- 
Fourier transform [75], whose inverse is in fact more complex 
than that of the standard Fourier transform (corresponding 
to q — 1). Indeed, it involves some subtleties presently under 
study. 

2 The (q, a)-distributions have, for q > 1 and a < 2, two 
distinct power-law regimes, namely the intermediate and the 



Putting all this together, we may alert the reader about 
an error commonly appearing in the literature. If (compu- 
tational, or experimental, or observational) data exhibit 
an asymptotic power-law behavior, many authors claim 
having shown the emergence of a Levy distribution. This 
is clearly incorrect: it could be a Levy distribution, or a 
q-Gaussian (which only concides with a Levy distribution 
for q = 2, in which case it corresponds to the a = 1 Levy 
distribution, i.e., the Cauchy-Lorentz distribution), or a 
(q, a)-stable distribution, or even none of them! Only the 
study of the central part of the empiric distribution (i.e., 
the non-asymptotic region) could in principle discriminate 
between all these possibilities. For example, Levy distri- 
butions with 1 < a < 2 always present an inflection point 
if represented in log-log scales; q-Gaussians do not have 
such inflection point in the same representation (see [731 
[74] for illustrations of this fact). 

The ubiquity of Gaussians in nature is generally be- 
lieved to be a consequence of the standard CLT. Analo- 
gously, one expects the ubiquity of q-Gaussian distribu- 
tions in asymptotically scale-invariant natural, artificial 
and social systems to be a consequence of the presently 
described q-generalized CLT. Two such systems are under 
active investigation nowadays, namely unimodal maps, 
and long-range classical Hamiltonians. 

Simple unimodal dissipative maps Xt+i = f{xt) (like 
the logistic and similar ones) exhibit, for control param- 
eter values such that the Lyapunov exponent is positive 
(i.e., strong chaos), sums of many successive iterates which 
approach Gaussians. If the control parameter is, instead, 
close to the edge of chaos (where the Lyapunov exponent 
is zero), q-Gaussians appear to be gradually approached 
by the sums of many successive iterates: see details in [521 

mm- 

Long-range-interacting many-body classical Hamilto- 
nians such as the HMF one [79] are known to present ther- 
mostatistical anomalies. Indeed, if isolated (i.e., micro- 
canonical ensemble), they can present, for a range of values 
for the total energy and some classes of initial conditions 
(e.g., the so-called water-bag and double- water-bag ones) 
a negative specific heat. Furthermore, they exhibit a max- 
imal Lyapunov exponent which vanishes in the N — > oo 
limit [54 , a fact which seemingly breaks down the validity 
of Boltzmann's molecular chaos hypothesis. In such a sys- 
tem, various other anomalies have been exhibited such as 
aging, spin-glass nature and non-Maxwellian distribution 
of velocities (see [80] and references therein). This state 
is usually referred in the literature as the QSS (quasi- 
stationary state) or metastable one. It has been recently 
exhibited that the system is nonergodic in the QSS, its 
ensemble-averaged and its time-averaged distribution of 
velocities being both non-Maxwellian and different among 
them. Its time-averaged one, which is the typically ob- 
served in experiments, appears to approach a q-Gaussian: 
see details in [T6 ] fT7 ] (see also [81(182]). 



distant ones (see details in [731174] ). The intermediate power- 
law approaches q-Gaussians in the a — » 2 limit. The distant 
one approaches Levy distributions in the q — > 1 limit. 
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Fig. 1. The function y = e~ x ^ with £ = 10 in typical repre- 
sentations: (a) linear-linear; (b) log-linear; (c) log-log. In the 
linear- linear representation, the two regimes x/£ — » and 
x/£ — > oo are clearly visible. Not so in the log-linear and log-log 
representations. 



3 Distinct physical mechanisms versus distinct 
regimes of a single physical mechanism 

We address in this section a question which emerges quite 
frequently in physics and elsewhere. It is in particular 
quite often present in the literature of experimental high 
energy physics. The question concerns within which con- 
text can/must we analyze the possible concurrence of var- 
ious distinct physical mechanisms in a single set of ex- 
perimental data or phenomenon. Let us first remind a 
historical illustration: electricity magnetism and optics 
were discussed during centuries as being different physi- 
cal mechanisms. But, since the celebrated unification per- 
formed in Maxwell equations, these phenomena are cur- 
rently thought as different asymptotic limits (or regimes) 
of a single mechanism, rather than as the coexistence 
of different physical mechanisms. The same happened of 
course in the case of the weak and electromagnetic in- 
teractions. Therefore, along such discussions we must al- 
ways have in mind that the opposition of different regimes 
against different physical mechanisms might, in this sense, 
be a historical and possibly provisory one. 

Let us mathematically illustrate this in another fa- 
mous example, namely Planck's law for the black-body 
radiation. It reads, for the photon energy density per unit 
volume, 

/ x 8-Khv 3 

= ^hv/ksT _ I) • ( 41 ) 

This law presents a low-frequency asymptotic regime, 

u{v) ~ 8irk B Tc- 3 v 2 , (42) 

as well as a high-frequency asymptotic regime 

u(v) ~ ^hc- 3 u z e- h ^ kBT , (43) 

with an intermediate or crossover regime in between (where 
u(v) is roughly constant). The high-frequency regime obvi- 
ously is a quantum one (indeed, we see therein Planck con- 
stant h), whereas the low- frequency regime is a classical 
one. During many years these two regimes were considered 
as two distinct physical mechanisms. Since the formulation 
of quantum mechanics, however, many scientists prefer to 
think of them in a more "modern" manner, namely as two 
distinct asymptotic limits of a single physical mechanism. 

In too many occasions, such discussions are regretfully 
driven, not by fundamental questions, but rather by the 
particular way of representing the results. For example, if 
we have the relation y = e~ x ^ with x > and £ > 0, we 
can of course distinguish two asymptotic regimes, namely 
the < x « £ and the x » £ ones. These two regions 
can be clearly identified in a linear-linear representation 
(y <~ 1 — (x/£) for x/t; — > 0, and y <~ for x/£ — ► oo), 
less clearly however in a log-linear representation (where 
we only see a single straight line), or in a log-log represen- 
tation (where we see a downwards ever bending curve): 
see Fig. 1. But, in such a situation, few scientists would 
keep talking of two distinct physical mechanisms; most of 
them would prefer to talk about two regimes of a single 
mechanism. 
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Fig. 2. The function y = e~ x/i = 1/[1 + (q- l)x/£] 1/(q - 1) with 
q = 2 and £ = 10 in typical representations: (a) linear-linear; 
(b) log-linear; (c) log-log. In the linear-linear, log-linear (see 
the departure from a straight line) and log-log representations 
(see the roughly horizontal and the power-law regions), the two 
regimes x/£ — > and x/£ — > oo are clearly visible. Not so in 
the (d) qlog-linear and (e) qlog-qlog representations, as can be 
seen in Fig. 3. 



Let us discuss now this very same point for a relation 
y = eq x/i = l/[l + ((7-l)x/';] 1 /(«- 1 ) withq > 1, which fre- 
quently emerges in high-energy collisions, and elsewhere. 
Once again, we can neatly distinguish two asymptotic 
regimes, namely the x/£ — > and the x/£ — > oo ones. 
In the first regime we have y ~ 1 — (x/£) (Vg > 1); in the 
second one, we have y ~ 0, or, in a more detailed form, 
y ~ l/[(q- ljar/f] 1 ^ 9-1 ) (Vg > 1). Let us now check pos- 
sible representations of this relation in Figs. 2 and 3. Like 
for the q = 1 previous example, we clearly see the two 
regimes in the linear-linear representation. The novelty 
for q > 1 is that we also see it in the log-linear and in 
the log-log representations. However, there is no a priori 
fundamental reason for now talking about two different 
mechanisms instead of just keep talking, as for the q = 1 
case, about two different asymptotic regimes. In fact, as 
can be seen in the qlog — linear representation, a single 
straight line represents the whole curve, in total similarity 
with the log — linear representation for the q = 1 case. 
Also, analogously to the q = 1 case, in a qlog — qlog rep- 
resentation we do not clearly see the two regimes. These 
few simple remarks might help to put in a proper context 
the discussion of experimental data, thus avoiding naive 
conclusions and enabling a deeper understanding of the 
physical phenomena which are involved. 



4 Conclusions 

The q-exponential functions, frequently present (at least 
as approximations) in high-energy physics and elsewhere, 
emerge naturally within nonextensive statistical mechan- 
ics by optimizing the nonadditive entropy S q under ap- 
propriate constraints. Their origin at a mesoscopic level 
might include multiplicative noise, non-Mar kovian pro- 
cesses, space-dependent diffusion coefficients, nonlinear dif- 
fusion, and related mechanisms. The discussion at this 
level will naturally provide further understanding but will 
be in general insufficient for determining the precise value 
of an index such as q. To achieve this deeper understand- 
ing, the analysis of a microscopic model (or, more pre- 
cisely, a class of models) will be necessary. Typically, if the 
microscopic dynamics exhibits, for classical systems, van- 
ishing Lyapunov exponents (or analogous situations for 
quantum systems), the system automatically becomes a 
strong candidate for the use of the nonextensive thermo- 
statistical concepts. Such situations might occur, for in- 
stance, whenever the generic correlations are long-ranged 
in space-time (long-ranged interactions, quantum strongly 
entangled systems, etc). By studying quantities such as 
the sensitivity to the initial conditions as a function of 
time, one may obtain, from first principles, the value of 
the index q (or of indices related to q) . 
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